RELATIVE OSCILLATION THEORY FOR DIRAC OPERATORS 
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Abstract. We develop relative oscillation theory for one-dimensional Dirac 
operators which, rather than measuring the spectrum of one single operator, 
measures the difference between the spectra of two different operators. This 
is done by replacing zeros of solutions of one operator by weighted zeros of 
Wronskians of solutions of two different operators. In particular, we show that 
a Sturm-type comparison theorem still holds in this situation and demonstrate 
how this can be used to investigate the number of eigenvalues in essential 
spectral gaps. Furthermore, the connection with Krein's spectral shift function 
is established. As an application we extend a result by K.M. Schmidt on the 
finiteness/infiniteness of the number of eigenvalues in essential spectral gaps 
of perturbed periodic Dirac operators. 



1. Introduction 

To set the stage, let I ~ {a,b) C R (with — oo < a < b < oo) be an arbitrary 
interval and consider the Dirac differential expression 

(1.1) T= i(72^+0(a;). 

1 ax 

Here 

(1.2) (j){x) = (pciix)! + (/>am(a;)CTi + {m + (l>sc{x))a3, 
(Ji, (72, denote the Pauli matrices 

(1.3) = ; J ) ' = ( i 7 ) , ^3 = ( J \ y 

and TO, (/)sc, ^ci, and i/iam are interpreted as mass, scalar potential, electrostatic 
potential, and anomalous magnetic moment, respectively (see [TH], Chapter 4). As 
usual we require to S [0, oo) and 0sc , 0ci , '/'am G ^lod^) real- valued. We don't 
include a magnetic moment f = r + a2(t>Tng{x) since it can be easily eliminated by 
a simple gauge transformation r = UfU~^, U = exp(i (l)^g{r)dr) (there is also 
a gauge transformation which gets rid of (jjam or c/fd (see [7], Section 7.1.1)). 

If r is limit point at both a and b, then t gives rise to a unique self-adjoint 
operator H when defined maximally (cf., e.g., [7], [H], [10] )■ Otherwise, we need 
to fix a boundary condition at each endpoint where t is limit circle. 

Explicitly, H is given by 

(1.4) H : D{H) L^{IX^) , 

/ ^ rf 
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where 

(1.5) D{H) = {fe L\I,C^)\ f e ACiocilX^), rf G L^I,C^), 

Wa{u-,f) = Wbiu+,f)=0} 

with 

(1.6) WAf.g) = i{nx),a2g{x)) = fi{x)g2{x) - h{x)gi[x) 

the usual Wronskian (we remark that the Hmit Wa,b{-, ■■) = limx-s-a.f) Wx{-, ■■) exists 
for functions as in (ll.5|) ). Here the function u- (resp. u+) used to generate the 
boundary condition at a (resp. b) can be chosen to be a nontrivial solution of 
Tu = if T is limit circle at a (resp. b) and zero else. 

We refer to the monographs [7| , [21] , [J^ for background and also for further 
information about Dirac operators and their applications. 

However, even though the Dirac operator is as important to relativistic quantum 
mechanics as the Schrodinger operator to nonrelativistic quantum mechanics, much 
less is known about its discrete spectrum. The main reason of course being that in 
contradistinction to typical Schrodinger operators, Dirac operators are not bounded 
from below and thus approaches relying on semi-boundedness are not applicable. 

Our aim in the present paper is to develop what we will call relative oscillation 
theory for a pair of Dirac operators Hi and associated with two potentials ipi and 
^0 as above. As we will show, it turns out to be an effective tool for both counting 
eigenvalues in essential spectral gaps as well as for investigation the accumulation 
of eigenvalues at the boundary of an essential spectral gap. 

Let {f,g) = flgi + /1 32 and |/| = + 1/2!^ denote the scalar product 

and norm in C^. Our key ingredient will be the Wronskian of two (nontrivial) 
real-valued solutions uq and ui satisfying tqUq = AqUq and tiUi = AiUi. Then we 
define a Priifer angle for the Wronskian W{ua, ui) via 

(1.7) ^1 = ^(-) ('"mS'^JO ■ 

^ ' \W,x{ui,-ia2Uo)J ^ \cos{i/j{x)) J 

Note that ■0(a;) is uniquely determined up to a multiple of 2tt by the requirement 
that ipix) should be continuous since the two Wronskians cannot vanish simulta- 
neously. 

The total difference 

(1.8) #(c.d) ("0, ui) = rAi,o(d)/7r] - LAi,o(c) AJ - 1 

will then be called the weighted number of sign flips of the Wronskian W(uo,ui) 
in the interval (c, d) <Z I (with a < c < d < b). Here \ x\ = max{n G Z | n < a;} and 
\x~\ = min{n £ Z\n > x} are the usual floor and ceiling functions. 

In fact, #(c,ci)(wo. Ml) counts the number of sign flips of W{uo,ui) where a sign 
flip is counted as -1-1 if i/; increases along the sign flip and as —1 if ■0 decreases. 
Moreover, one can show that a zero xq is counted as +1 if {uq{xq), A(j){xQ)uo{xo)) > 
and as —1 if (mo(xo), A(/)(xo)uo(xo)) < 0, where 

(1.9) A0 -01-00. 
We will also set 

(1.10) ^{uo,ui)^ lim 4ic4){uq,ui) 

provided this limit exists. This will for example be the case if the perturbation is 
of a definite sign, A0(a;) > or A0(a;) < 0, at least for x near a and b. We will call 
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Ti — Ai relatively nonoscillatory with respect to tq — Xq if #(uo,ui) is finite and 
relatively oscillatory otherwise. 

Our first result implies that if we choose uq and ui to be Weyl solutions, then the 
weighted number of sign flips counts precisely the eigenvalue difference. Recall that 
a solution m_(z, .) of tu — zu is called Weyl solution at a if it is square integrable 
near a and fulfills the boundary condition of H at a (if there is any, i.e., if r is 
limit circle at a). Such a solution is unique up to a constant if it exists (e.g. if 
z ^ aess{H)) and it can be chosen to be real for z G M. Similarly a Weyl solution 
M+(z, .) at b is defined. 

Finally, denote by Pn{H), C R, the family of spectral projections associated 
with the self-adjoint operator H (see e.g. [T8j). 

Theorem 1.1. Let Hq, Hi be self-adjoint operators associated with tq, ti, respec- 
tively, and separated boundary conditions. Suppose 

(i) A</) < 0, near singular endpoints, 

(ii) liuix^a A0(a;) ~ if a is singular and Imix^b A(f)(x) = if b is singular, 

(iii) Hq and Hi are associated with the same boundary conditions near a and 
b, that is, Mo,- (A) satisfies the boundary condition of Hi at a (if any) and 
Mi,+ (A) satisfies the boundary condition of Hq at b (if any). 

Suppose Xq < inf (Jess (i/o). Then 

(1.11) dimRanP(_oo,Ao)(-H^i) - dimRanP(_oo.Ao](Jyo) = #(wi,t(^o), wo,±(Ao)). 
Suppose aess{Ho) H [Ao,Ai] = 0. Then ti — Aq is relatively nonoscillatory with 

respect to tq — Ao and 

dim Ran P[Ao , Ai ) [Hi ) - dim Ran P( Ao , Ai] [Hq ) 

(1.12) =#(ui,t(Ai),uo,±(Ai)) -#(mi,^(Ao),mo,±(Ao)). 

The proof will be given at the end of Section (5] 

Remark 1.2. Note that condition (ii) implies aessiHo) = o'essiHi) (cf. Lemma \2. 7| 
below). In addition, (ii) implies that any function which is in ©(tq) near a (or b) 
is also in S)(ti) near a (or b), and vice versa. Hence condition (iii) is well-posed. 

In the case where the resolvent difference of Hi and Hq is trace class, the differ- 
ence in p.l2p as opposed to p. lip can be avoided if we replace the left-hand side by 
Krein's spectral shift function ^(A, Hi, Hq) (see [23] for more information on Krein's 
spectral shift function). In order to fix the unknown constant in the spectral shift 
function, we will require that Hq and Hi are connected via a path within the set of 
operators whose resolvent difference with Hq are trace class. Hence we will require 

Hypothesis 1.3. Suppose Hq and Hi are self-adjoint operators associated with tq 
and Ti and separated boundary conditions. Assume that 

• A(f) is relatively bounded with respect to Hq with HQ-bound less than one, 
and 

• y^\A(j)\{HQ — z)^^ is Hilbert-Schmidt for one (and hence for all) z e p{Hq). 

It was shown in |B1 Sect. 8] that these conditions ensure that we can interpolate 
between Hq and Hi using operators 7?^, e G [0, 1], such that the resolvent difference 
of Hq and H^ is continuous in e with respect to the trace norm. Hence we can fix 
^(A, Hi,Hq) by requiring e >-)• ^(A, H^,Hq) to be continuous in L^(R, (A^ -I- l)~^dX), 
where we of course set ^(A, Hq, Hq) = 0. While ^ is only defined a.e., it is constant 
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on the intersection of the resolvent sets R n p{Ho) p{Hi), and we wiU require it 
to be continuous there. In particular, note that by Weyl's theorem the essential 
spectra of Hq and Hi are equal, aess{Ha) = <^ess{Hi). Then we have the following 
result: 

Theorem 1.4. Let Hq, Hi satisfy Hypothesis ] 1.3[ Then for every A G Rrip{Ho)r\ 
p{Hi) we have 

(1-13) a^,HuHo) = #(V-'o,±(A),Vi,^(A)). 

Again, the proof will be given at the end of Section [5] 

In particular, this result implies that under these assumptions ti — A is relatively 
nonoscillatory with respect to tq — A for every A in an essential spectral gap. 

Concerning the history of these results we mention that the analogs of Theo- 
rem 11.11 and Theorem 11.41 were first given in the case of Sturm-Liouville operators 
by Kriiger and Teschl [6] , [4] extending earlier work of Gesztesy, Simon, and Teschl 
[3] which corresponded to the case Hi = Hq. In the case of Dirac operators the 
case Hi = Hq was first given in Teschl ^7] . 

Finally, we will show how #(mo,wi) can be used to settle the question whether 
the eigenvalues introduced by a given perturbation will accumulate at a boundary 
point of the essential spectrum and apply this to the case of perturbed periodic 
Dirac operators. 

We first recall some basic facts from the theory of periodic Dirac operators (cf., 
e.g., [21], Chapter 12, [21], Chapter 16). Let Hq be a Dirac operator associated with 
periodic potential (po of period a > 0, that is, (fioix + a) = </>o(a;), x G I = {a, oo). 
The essential spectrum of Hq is purely absolutely continuous and consists of a 
countable number of bands, that is, 

(1.14) aessiHQ)=\J[E2j,E2j + l] 

jez 

with • • • E2j < E2j+i < E2j+2 < E2j+3 ■ ■ ■ ■ In addition, in every essential spectral 
gap there can be at most one eigenvalue. 

Moreover, Floquet theory implies the existence of an (anti-)periodic solution 
UQ{Ej,x) at each boundary point of the essential spectrum. 

To phrase our result, we recall the iterated logarithm log„(a;) which is defined 
recursively via 

logo(a;) = X, log„(a;) = log(log„_i(a;)). 
Here we use the convention log(a;) = log \x\ for negative values of x. Then log„(x) 
will be continuous for x > e„_i and positive for x > e„, where C-i = — cxd and 
Cra = e°"^^ Abbreviate further 

1 " 

Explicitly we have 

Lq{x) ^ X, Li{x) ^ x\og{x), L2{x) ^ x\og{x)\og{\og{x)), ... 
With this notation we have the following result: 

Theorem 1.5. Let Ej be a boundary point of the essential spectrum of the periodic 
operator Hq and let uq(x) be a corresponding (anti-) periodic solution of tqUq = 
EjUQ. 



RELATIVE OSCILLATION THEORY FOR DIRAC OPERATORS 



5 



Suppose 

n 

(1.15) M^) = M^) - 4 E rr^*^^'" + o(L„(x)-2) 

for some constant matrices (f>i,k, < k < n, and define 

- _ 2_ {u{x), ((m + (f)o.sc{x))cr3 + (l)o.i^m{x)(Ti)u{x)) 

(1.16) Bk^ / {u{x),(l)i ku{x))dx, 0<k<n. 

a Ja 

Then the eigenvalues of Hi accumulate at Ej if 

(1.17) ABa = • • • = ABn-i = 1 and AB^ > 1 
and the do not accumulate at Ej if 

(1.18) ABq = ■■■ = ABn-i = 1 and AB„ < 1. 
The proof will be given at the end of Section |4l 

In the case of Sturm-Liouville operators this result originates in the work of 
Rofe-Beketov [H]-[II] (see also the recent monograph ^3j) who proved the case 
n — 0. His work was recently improved by Schmidt [15] who gave a new proof and 
obtained the cases n = 0, 1. Extending the approach by Schmidt the general case 
was obtained in Kriiger and Teschl [5^. Schmidt also established the case n = 0,1 
for Dirac operators in [16] . In his paper [16] he also gives an equivalent formulation 
for the criterion in terms of the gradient of the Floquet discriminant and shows how 
the above criterion can be applied to radial Dirac operators via a transformation 
from [14]. In fact, if 

1 d k 

(1-19) Tk^-a2— + -a3 + (t)ir), re(0,cx)), 

1 ar r 

is a radial Dirac operator (i.e. one which arises by separation of variables in spherical 
coordinates \19i Sect. 4.6.6]), then the unitary transformation (^14j Lem. 3]) 

transforms t to 

(1-21) U*tU ^la,i-+\\ll+'--l]a, + ^p^l + <t>{r). 

Since 

(1-22) [f^ - l) -3 + ^(^1 = + 1)^ + Oir-^) 

our result is directly applicable to this situation. 

We also refer to [16] and the recent work by Cojuhari [2] for more on the history 
of this problem and references to related results. Analogous results for the discrete 
case, Jacobi matrices, can be found in [1]. 
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2. Relative Oscillation Theory 

After these preparations we are now ready to develop relative oscillation theory. 
Our presentation closely follows [S]- 

Definition 2.1. For tq, ti possibly singular Dirac operators as in il.l]) on (a, 6), 

we define 

(2.1) #(uo,Mi) = liminf #(c,rf)(Mo,'Ui) and #(wo, "i) = limsup #(c,rf)(wo, "i), 

where TjUj — XjUj , j = 0,1. 

We say that ^{uq,ui) exists, if ^{uo,ui) = #(mo,wi), and write 



in this case. 

By Lemma 13.11 below one infers that #(uo,mi) exists if (j)o ^ Xo ^ + Xi has 
the same definite sign near the endpoints a and b. On the other hand, note that 
#(uo,ui) might not exist even if both a and b are regular, since the difference of 
Priifer angles might oscillate around a multiple of tt near an endpoint. Furthermore, 
even if it exists, one has #(wo, ui) ~ if(a,b)iuo, ui) only if there are no zeros at the 
endpoints (or if — Aq — + Ai > at least near the endpoints). 

We begin with our analog of Sturm's comparison theorem for zeros of Wron- 
skians. We will also establish a triangle-type inequality which will help us to provide 
streamlined proofs below. Both results follow as in fB]. 

Theorem 2.2 (Comparison theorem for Wronskians). Suppose Uj satisfies TjUj = 
XjUj, j ~ 0, 1, 2, where Aq — 0o < Ai — (/>i < A2 — 02- 

If c < d are two zeros ofWx{uo,ui) such that Wx{uo,ui) does not vanish identi- 
cally, then there is at least one sign flip of Wxiuo,U2) in (c, d). Similarly, if c < d 
are two zeros of Wx{ui,U2) such that Wx{ui,U2) does not vanish identically, then 
there is at least one sign flip of Wx{uq,U2) in (c, d). 

Theorem 2.3 (Triangle inequality for Wronskians). Suppose Uj, j — 0,1,2 are 
given real-valued non-vanishing vector functions. Then 

(2.3) #(uo,Ui) +#(U1,U2) - 1 <±{UQ,U2) < #(mo,wi) +#(mi,?^2) + 1, 

and similarly for f^ replaced by f/=. 

Definition 2.4. We call ti relatively nonoscillatory with respect to tq, if the quan- 
tities ^{uq, Ui) and #(wo. Mi) are finite for all solutions tjUj = 0, j = 0, 1. 

We call Ti relatively oscillatory with respect to tq, if one of the quantities #(wo, Wi) 
or ^(wo,Mi) is infinite for some solutions tjUj = 0, j — 0, 1. 

Note that this definition is in fact independent of the solutions chosen as a 
straightforward application of our triangle inequality (cf. Theorem 12. 3p shows. 

Corollary 2.5. Let TjUj — TjVj = 0, j = 0, 1. Then 



The bounds can be improved using our comparison theorem for Wronskians to 
be < 2 in the case of perturbations of definite sign. 

To demonstrate the usefulness of Definition 12.41 we now establish its connection 
with the spectra of self-adjoint operators associated with t,-, j — 0, 1. 



(2.2) 



#{uo,ui) = #(ito,ui) = #(-uo,wi) 



(2.4) 



\#{uo,ui) - #{vo,vi)\ <4, |#(wo,wi) -#(wo,«i)| <4. 
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Theorem 2.6. Let Hj be self-adjoint Dirac operators associated with Tj, j = 0, 1. 
Then 

(i) To — Ao is relatively nonoscillatory with respect to tq — Ai if and only if 

dimRanP(Ao,Ai)(-ffo) < oo. 

(ii) Suppose dimlian P(^Xg x^^{Hq) < oo and ti ~ X is relatively nonoscillatory 
with respect to tq — X for one X G [Aq, Ai]. Then it is relatively nonoscilla- 
tory for all A e [Ao,Ai] if and only if dimKan P^Xo.\i)iHi) < 

Proof. Item (i) is [T7J Thm. 4.5] and item (ii) follows as in ^64. □ 

For a practical application of this theorem one needs of course criteria when 
Ti — A is relatively nonoscillatory with respect to tq — A for A inside an essential 
spectral gap. 

Lemma 2.7. Let limx^a{4>o{x) ~ = if a is singular, and similarly, 

lim.x^bi'Poix) — 01 (2;)) =0 if b is singular. Then aessiHo) = o'essiHi) and ti — X 
is relatively nonoscillatory with respect to tq — X for X € M\aessiHo)- 

Proof. Since ti can be written as ti = tq + c/jq + cjji, where (j)Q has compact support 
near singular endpoints and |0i| < e, for arbitrarily small e > 0, we infer that 
Rhi{z) — Rh„{z) is the norm limit of compact operators. Thus Rhi{z) — Rhq{z) 
is compact and hence aess{Ho) = aess{Hi). 

Let i5 > be the distance of A to the essential spectrum and choose a < c < d < 6, 
such that 

\cj)i{x)-cj)o{x)\<5/2, x(^{c,d). 
Clearly #(c,rf)('UOj "^ti) < 00, since both operators are regular on (c, d). Moreover, 
observe that 

0O-A+ <0i-A<0o-A_, A±=A±<5/2, 
on / = (a, c) or / = (d, b). Then Theorem 12.61 (i) implies #/(uo(A-), uq{X+)) < 00 
and invoking Theorem 12.21 shows #/(uo(A±), mi(A)) < 00. From Theorem 12.31 and 
12.61 (i) we infer 

#^(uo(A),wi(A)) < #/(uo(A),uo(A+)) + #/(uo(A+), ui(A)) + 1< 00, 

and similarly for #^(uo(A), ui (A)). This shows that ti— A is relatively nonoscillatory 
with respect to tq. □ 

Our next task is to reveal the precise relation between the number of weighted 
sign flips and the spectra of Hi and Hq. The special case Hq = Hi is covered by 

Theorem 2.8 f |17i Thm. 4.5]). Let Hq be a self-adjoint operator associated with 
To and suppose [Aq, Ai] fl aessiHa) = 0. Then 

(2.5) dimRanP(A„^A,)(i/o) = #(V'o,t(Ao), V^o,±(Ai)). 

Combining this result with our triangle inequality already gives some rough 
estimates in the spirit of Weidmann [5D] who treats the case Hq — Hi . 

Lemma 2.9. For j = 0,1 let Hj be a self-adjoint operator associated with Tj and 
separated boundary conditions. Suppose that (Ao,Ai) C M\{aess{Ho) U aessiHi)) , 
then 

dimRanP(Ao,Ai)(^i) - dimRanP(Ao,Ai)(^o) 

(2.6) <#(^i.t(^i)'V'o.±(Ai))-#(Vi,t(^o),V'o,±(Ao)) + 2, 
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respectively, 

diniRanP(Ao,Ai)(^i) - dimRanP(;^Q_;^^)(i7o) 

(2.7) > #(^i,t(Ai), V^o.±(Ai)) - #(V^i,t(Ao), Vo,±(Ao)) - 2. 

Given these preparations the proofs of Theorem 11.11 and Theorem 11.41 ean be 
done as in [6]. 

Proof of Theorem For the proof one can hteraUy follow the arguments in Sec- 
tion 6 of [B]. The only noteworthy difference is that in Lemma 6.4 one has to use the 
limsup of the largest eigenvalue and the liminf of the lowest eigenvalue of ^. □ 

Proof of Theorem \1.4\ For the proof one can literally follow the arguments in Sec- 
tion 7 of 0. □ 

3. More on Prufer angles and the case of regular operators 

The purpose of this section is to collect some further facts on Priifer angles for 
Wronskians and to prove Theorem 1 1.1 1 in the case of regular operators. Even tough 
the Priifer angle Ai^o introduced below is different from ip used in the introduction 
it will be equivalent for our purpose (cf. Definition 14.11 below) . We closely follow 
[B] and we will provide proofs only when there is a significant difference to the 
Sturm-Liouville case. 

We first introduce Priifer variables for u G C(/,R^) defined by 

(3.1) ui{x) ^ pu{x) sin{Ou{x)) U2(x) = cos(6'„(x)). 

If u is never (0,0) and u is continuous, then p„ is positive and 9u is uniquely 
determined once a value of Ou{xo), Xq & I is chosen by the requirement 9u G C{I, R). 
The connection with the Wronskian is given by 

(3.2) Wa;(M,w) = -p„(x)pi,(x)sin(A„,„(x)), A^,„(a;) = 6'^(x) - 6'„(x). 

Hence the Wronskian vanishes if and only if the two Priifer angles differ by a 
multiple of tt. We will call the total difference 

(3.3) [uo, ui) = rAi,o(d)/7rl - LAi,o(c) AJ - 1 

the number of weighted sign flips in (c, d), where we have written Ai,o(a;) = Au^^uo 
for brevity. 

Next, let us take two real- valued (nontrivial) solutions Uj, j = 1, 2, of TjUj = XjUj 
and associated Priifer variables pj , 9j . Since we can replace — A it is no 

restriction to assume Ao = Ai = 0. 

Under these assumptions Wx{uo,ui) is absolutely continuous and satisfies 

(3.4) W!^{uo,ui) = {uo{x), {(j)i{x) - <j)oix))ui{x)) . 

Lemma 3.1. Abbreviate Ai_o(x) = 0i{x) — 0o{x) and suppose Ai_o(a;o) = mod tt. 
If —{uo{x), A(p{x)ui{x)) is (i) negative, (ii) zero, or (Hi) positive for a.e. x G 
(xo,a;o -f e) respectively for a.e. x G (xq — £tXq) for some e > 0, then the same is 
true for (Ai_o(a;) - Ai^o(a;o))/(2; - a;o). 

Hence i^[c.d)i'u-o,ui) counts the weighted sign flips of the Wronskian Wxiuo,ui), 
where a sign flip is counted as +1 if — A0 is positive in a neighborhood of the 
sign flip, it is counted as —1 if — A0 is negative in a neighborhood of the sign flip. 
If A(j) changes sign (i.e., it is positive on one side and negative on the other) the 
Wronskian will not change its sign. In particular, we obtain: 
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Lemma 3.2. Letuo, ui solve TjUj = 0, j = 0, 1, where Aip < 0. Then i^(^a,b)(,uo,ui) 
equals the number sign flips of W{uo,ui) inside the interval (a, 6). 

In the case A(j) > we get of course the corresponding negative number except for 
the fact that zeros at the boundary points arc counted as well since [— .xj = — [x]. 
That is, if A0 < 0, then #(c,d)(Mo, ^i) equals the number of zeros of the Wronskian 
in (c, d) while if > 0, it equals minus the number of zeros in [c, d]. In the next 
theorem we will see that this is quite natural. In addition, note that jf{u, u) = — 1. 

Finally, we establish the connection with the spectrum of regular operators. A 
finite end point is called regular if all entries of (j) are integrable near this end 
point. In this case boundary values for all functions exist at this end point. In 
particular, r is called regular if both end points a, b are regular. In the regular case 
the resolvent of H is Hilbert-Schmidt and hence the spectrum is purely discrete 
(i.e., aess{H) = 0). 

Theorem 3.3. Let H^, Hi be regular Sturm Liouville operators associated with 

To, Ti and the same boundary conditions at a and b. Then 

(3.5) 

dim Ran P(_oo,Ai)(-ffi) - dim Ran P(_cx>,Ao](-f^o) = #(o,6)(wo,±(Ao), wi,t(Ai)). 

The proof will be given below employing interpolation between Hq and Hi, using 
Hs = {1 — e)Ho + eHi together with a careful analysis of Priifer angles. 

It is important to observe that in the special case Hi ~ Hq, the left-hand side 
equals dimRan P(^Xi,Xo){^o) if Ai > Aq and — dimRan P[x„^\j^]{Ho) if Ai < Aq. This 
is of course in accordance with our previous observation that #(wo,±(Ao), wi,:p(Ai)) 
equals the number of zeros in (a, b) if Ai > Aq while it equals minus the numbers 
of zeros in [a, b] if Ai < Aq. 

Now let us suppose that to,i are both regular at a and b with boundary conditions 

(3.6) cos(Q)/i(a)-sin(a)/2(a)=0, cos(/3)/i(6) - sin(;3)/2(6) = 0. 

Hence we can choose u±{\,x) such that w_(A, a) = (sin(a), cos(a)) respectively 
M+(A, 6) = (sin(/3), cos(/3)). In particular, wc may choose 

(3.7) 6'_(A,a) = a € [0,7r), -6'+(A, &) = tt - /3 € [0, tt). 
Next we introduce 

(3.8) Te = To + e{<j)i - <j)o) 

and investigate the dependence with respect to £ S [0,1]. If solves T^Ue = 0, 
then the corresponding Priifer angles satisfy 

(3.) ««)-^. 
where the dot denotes a derivative with respect to e. 
Lemma 3.4. We have 

(3.10) w.i-.±^u,^)={ ^^Pr^'Y^^^^^^ . 

t -Ja {ue,-{r),{Mr) - Mr))us,-{r))dr 
where the dot denotes a derivative with respect to s and Ue,±{x) = U£^±{0,x). 
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Denoting the Priifer angles of U6,±(x) — Ue.±{0jx) by 9e.+{x), this result implies 
for 00 - 01 > 0, 

de,+ {x) = ^ < 0, 



(3.11) e,A^)^ '-' ' "''"''^ ' ^j;' " ' " >o, 



XfK.-W, ( Mr) - (t)i {r))ue,-{r))dr 

with strict inequalities if 4>q > 0i on a subset of positive Lebesgue measure of (x, 6), 
respectively (a, x). 

Now we are ready to investigate the associated operators Hq and Hi . In addition, 
we will choose the same boundary conditions for as for Ho and Hi. 

Lemma 3.5. Suppose (/)o~0i ^0 (resp. 00 ~ 0i < OJ. Then the eigenvalues of H^ 
are analytic Junctions with respect to e and they are decreasing (resp. increasing). 

In particular, this implies that dimRanP(_oo,A)(-ffe) is continuous from below 
(resp. above) in e if 0o — 0i ^0 (resp. 0o — 0i < 0). 
Now the proof of Theorem 13.31 follows literally as in [5] . 

4. Relative oscillation criteria 

As in the previous sections, we will consider two Dirac operators tj, j = 0,1, 
and corresponding self-adjoint operators Hj, j = 0,1. Now we want to answer 
the question, when a boundary point E of the essential spectrum of Hq is an 
accumulation point of eigenvalues of Hi. By Theorem 12.61 we need to investigate if 
Ti — E is relatively oscillatory with respect to To — i? or not, that is, if the difference 
of Priifer angels Ai o = 6i — 6q is bounded or not. 

Hence the first step is to derive an ordinary differential equation for Ai^o- While 
this can easily be done by subtracting the differential equations for 6i and 0Oj the 
result turns out to be not very effective for our purpose. However, since the number 
of weighted sign flips #(c,d)(Mo, wi) is all we are eventually interested in, any other 
Priifer angle which gives the same result will be as good: 

Definition 4.1. We will call a continuous function -0 a Priifer angle for the Wron- 
skian W{uo,ui), «/ #(c,<i)(Mo, ^i) = \ip{d)/Tr'] - [ip{c)/Tr\ - 1 for any c,de {a,b). 

Hence we will try to find a more effective Priifer angle ip than Ai^o for the 
Wronskian of two solutions. The right choice for Sturm-Liouville equations was 
found by Rofe-Beketov [8] (see also the recent monograph [13]) and it turns out the 
analogous definition is also the right one for Dirac operators |16) : 

Let uq, vq be two linearly independent solutions of {tq — X)u = with W{uq, vq) = 
1 and let ui be a solution of (ti — X)u — 0. Define "0 via 

(4.1) W{uo,ui) ^ -Rsm{^/j), W{vo,ui) ^ -Rcosiip). 

Since W{uo,ui) and W{vo,ui) cannot vanish simultaneously, -0 is a well-defined 
absolutely continuous function, once one value at some point xq is fixed. 



Lemma 4.2. The function ip defined in (|4.1[) is a Priifer angle for the Wronskian 
W{uo,ui). 

Proof. Since W{uo,ui) = —Rsm{ip) — —puoPui sin(Ai_o) it suffices to show that 
■0 — Ai^o mod 27r at each zero of the Wronskian. Since we can assume Ovq — 
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9uo G (0,7r) (by W{uq,vo) — 1), this follows by comparing signs of Rcos{^p) — 
)• □ 

Lemma 4.3. Let uq, vq be two linearly independent solutions of (tq — X)u = with 
W{uo,vo) = 1 and let ui be a solution of (ri — X)u = 0. 

Then the Priifer angle ip for the Wronskian W{uq, ui) defined in (I4.ip obeys the 
differential equation 

(4.2) ip' = -(uo cos(^) - Wo sin(V'), A(/){uo cos(-0) - vq sin(-0))), 

where 

/^4> = 4>i- (/)o- 

Proof. Observe R'ip' = -M^(uo, mi)' cos(V') + H/"(uo, wi)' sin(V') and use (15^ . (liTT|) 
to evaluate the right-hand side. □ 

To proceed we will need the following formula for a second solution of a Dirac 
equation which can be verified by a straightforward calculation: 

Lemma 4.4 fll2|. |16[ Lem. 1]). Let u be a nontrivial solution of tu — zu and 
choose xq E L Then 

I"' {u(r ),4:{r)u{r)) _ ■ 

|u(r)|4 '^'^ \u{x) 

where 

(4.4) (t>{x) = (to + 4>sc{x))a'i + (t>a.Ya{x)ai, 

is a second linearly independent solution satisfying W{u^v) = 1. 

Now we will choose vq to be given by (|4.3p and, following Schmidt |16j . perform 
a Kepler transformation 

(4.5) cot(^(x)) = i f cot(^(.)) - 2 r iMr)Mr^uo{r)) \ 



(4-3) v{x)^[2i 



Hr)V 

to obtain 

i( N 1 /r,('^o(a;),0o(a;)uo(a;)) ■ 2, r v-, , ■ , , v-, 1 1 w 
ip {x) =- 2 ^ — — sm {if{x)) + sin((^(a;)) cos{ip{x))- 

((cos(<^(.))-i^-^a2juo(x), 

(4.6) A0(x) ( cos(v5(x)) - -^-p^a^ u^{x)) 

Here we assume that a > is regular and 6 = cxd without loss of generality. Un- 
der the further assumption that |Mo(a;)|, |uo(x)|~^, and x^A0(x) are bounded this 
simplifies to 
(4.7) 

<pj{x) = - {a{x) sv[?(ip(x)) + sin(.^(a;)) cos(.^(a;)) -f B(x) cos^((p(x))) -t- 0(a;"^), 
where 

(4.8) A{x)=2 ■ — and B{x) = -{uo{x),x^ i:^(i)[x)uo(x)) . 
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Now we turn to the case where (jjoix) is periodic with period a > and choose uo 
to be the (anti-)periodic solution at a band edge. Taking averages 

(4.9) Jp(x) = — / ip(r)dr 

the above differential equation turns into (see [5] Section 5]) 

if- 

X 

where 



(4.10) Tp'{x) = ^(^Asin^(<^(a;))+sin(^(x))cos(^(x)) + B(a;)cos^(^(a;))^+0(a;"^), 



— ^ 2_ {uQ{x),(j)Q{x)uo{x)) 
a Jo l"o(a:^)|^ 



1 



X^OL 



(4.11) B{x) = --i (iio(r),r2A0(r)iio(r)). 
Moreover, if (/)i(x) is given by (|1.15l) then one computes 

1 " 2 

(4.12) B{x) = -4 E lT^^*^ + o{xH„{x)-'). 
Now we use the following result: 

Lemma 4.5 (|5l Lemma 4.7]). Fix some n £ No, let Q be locally integrahle on 
(a, 00) and abbreviate 



1 ""^ 1 



2 ■ 



T/ien solutions of the differential equation 
(4.13) ^ 

i^'(x) = - (sin^((p(x)) + sin(93(x)) cos((^(x)) - x^Q{x)cos^(if{x))) + of . „ ) 
fend to cxd if 



limsupL„(a;)^ {Q{x) - Q„(a;)) < -i 

a:— f 00 4 

and are bounded from above if 

liminf Ln{xf {Q{x) - Qn{x)) > -7. 

X^OO 4 

/n the last case all solutions are bounded under the additional assumption Q = 
g„(x)+0(L„(x)-2). 

Now this lemma implies Theorem 11.51 if A ~ 1. However, if A > we can easily 
reduce it to the case A = 1 by the simple scaling uo{x) — > (A)^/^ito(a;). which 
renders A — 1 and Bk — t- ABk- Similarly, if < we can reduce it to the case 
A > via the transformation (p —ip which renders A —A, Bk —Bk- Finally, 
in the case A = the result follows by using Proposition 1 from [16] (Lemma 5.1 
in [5]) in place of the above lemma. 
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